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Abstract. The correlation between baryon number and strangeness is used 
^ ' to discern the nature of the deconfined matter produced at vanishing chemical 

O potential in high-energy nuclear collisions at the BNL RHIC. Comparisons of 
results of various phenomenological models with correlations extracted from 
lattice QCD calculations suggest that a quasi-particle picture applies. At finite 
baryon densities, such as those encountered at the CERN SPS, it is demonstrated 
that the presence of a first-order phase transition and the accompanying 
development of spinodal decomposition would significantly enhance the number 
" ^ , of strangeness carriers and the associated fiuctuations. 
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l/^ , 1. Introduction 

o 

The goal of high-energy heavy-ion collisions is the creation of a new state of strongly 
interacting matter at very high energy density where the degrees of freedom carry 
color charges. These may be quark and gluon quasi-particles colored excitations 
of hadrons, or more complicated composite structures. The creation of such a state 
of matter is ushered in by the high energy densities created in such collisions. The 
reconversion of such an excited state of matter back into ordinary hadrons necessarily 
requires the presense of a phase transformation. The presence of such an excited state 
. of matter and a phase change between it and ordinary matter has been observed in 

' lattice simulations of QCD (LQCD) at a temperature of 170 MeV and vanishing baryon 

density [2]. The suppression of the number of high-momentum particles originating in 
the fragmentation of hard partonic jets requires the presence of a near opaque medium 
with densities at least 30 times that of conventional nuclear matter P]. Such densities, 
when compared with the results from LQCD require the system to lie firmly in the 
deconfined state. 

The phase transformation at a vanishing baryon density has been ascertained to 
be continuous. As a result, one would expect that the various conserved quantities 
such as the baryon number, charge and strangeness remain more or less unchanged 
within various sub-volumes of the phase space occupied by the system. As a result 
event-by-event fluctuations of these quantities may be used as probes of a state prior 
to hadronization. This will be the topic of the first part of these proceedings. It will be 
demonstrated that different assumptions regarding the nature of the Baryon number 
and strangeness carrying degrees of freedom lead to different measurable observables. 
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At zero temperature, most models predict the occurrence of a first-order phase 
transition when the density is raised. To date, there exist no lattice calculations in 
this regime. Thus, if there exists a first-order phase transition at zero temperature 
and finite chemical potential, this transition line should extend into the region of finite 
temperature and would persist until the chemical potential has dropped below a certain 
critical value /Xc 0]. There exist lattice QCD estimations |S] that suggest the presence 
of a first-order phase transition line at finite chemical potential as well as the existence 
of a critical end-point, though its precise location is not well determined. The second 
part of these proceedings will deal with how certain observables such as the ratio of 
strange to non-strange particles and the fluctuations of such ratios are influenced by 
the non-monotonic behaviour of the thermodynamic potential as the system traverses 
a first-order phase transition line in the course of its dynamical expansion. 

2. Baryonic strangeness at RHIC: /i << T 

Prior to the appearance of data from RHIC, it was believed that results from LQCD 
could eventually be explained within a picture of the deconfined phase as a weakly 
interacting plasma of quasiparticles The large radial and elliptic flow of the bulk 
matter however is somewhat inconsistent with the picture of a weakly interacting 
plasma [7|. Such observables seem to require the presense of a strongly interacting 
liquid. Indeed, recent results from lattice QCD calculation on spectral functions 
suggest the presence of bound, color-neutral states above Tc- This has led to the 
suggestion that at moderate temperatures, T ~ 1 — 2Tc, the system is composed 
of medium-modified (massive) quarks and gluons together with their (many) bound 
states [S]. In yet another picture, the system may be composed of strings between 
colored particles and their topological excitations [S| . At the present time, there exists 
no clear concensus on the nature of the degrees of freedom of the excited state of matter 
being created in RHIC collisions. 

Imagine that highly excited matter has indeed been created in the mid rapidity 
region at RHIC. The matter is strongly interacting and as a result has thermalized in 
the course of its expansion. The thermodynamics of this system may, at this point, 
be described in terms of a temperature T and chemical potentials iib, IJ-St IJ-Q for 
each of the conserved quantities of baryon number, strangeness and electric charge 
|1(J|. Experimentally determined chemical frcezeout temperatures place the chemical 
freeze out of the expanding gas of hadrons very close to the expected phase transition 
temperature. These features lead to the picture that the lifetime of the system as a 
gas of hadrons where the interactions are strong enough to change the chemistry of the 
state is rather short. It is this possibility which forms the basis of the assertion that 
the conserved quantities of baryon number B, strangeness S and electric charge Q, 
within a wide rapidity bin will receive minor contamination from neighbouring bins. 

As a result, the B,S and Q established in a given rapidity bin in the plasma 
phase will be approximately maintained by the hadron gas phase in each event. It 
may be argued that the correlation between the strangeness S and the baryon number 
B provides a useful diagnostic for the presence of strong correlations between quarks 
and anti-quarks. Considering a situation in which the basic degrees of freedom are 
non- interacting quarks and gluons, the strangeness is carried exclusively by the s 
and s quarks which in turn carry baryon number in proportion to their strangeness, 
Bg = —■^Sg, thus rendering strangeness and baryon number strongly correlated. This 
feature is in stark contrast to a hadron gas in which the relation between strangeness 
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and baryon number is less intimate. For example, at small chemical potential and 
temperature, the strangeness is carried primarily by kaons which have no baryon 
number. 

These considerations lead us to introduce the correlation coefficient |14) . 

_ , {BS)-{B) jS) (BS) 

^BS = ^•J 9 — — — o ,^„, . (il 

The last expression uses the fact that (5) vanishes. We have chosen the normalization 
such that Cbs is unity when the active degrees of freedom are individual quarks. 

When the system consists of independent species k having baryon number B^. and 
strangeness Sk, its total baryon number is -B = '^^nkBk and its total strangeness 
is S' = X^^rifcS'fc. The correlation coefficient may then be expressed in terms of the 
multiplicity variances af. = (nQ — {nkf^ w (n^), 

r< _ o Ek<^lBkSk _ o Efc {nk)BkSk 

~ T.k-isi ^ EkMsi ■ 

Thus, in the hadronic gas, the numerator receives contributions from only (strange) 
baryons (and anti-baryons), while the denominator receives contributions also from 
(strange) mesons, 

^ , (A) + (A) + ... +3(fl-)+3(n+) 

~ KO+KO + ... +9{n-)+9{Q+) ■ 

The hadronic freeze-out value value of Cbs is shown in Fig. ^ At the relatively 
high temperatures relevant at RHIC, the strange mesons significantly outnumber the 
strange baryons, so Cbs is smaller than unity. Indeed, including hadrons up to r2~, 
we find Cbs = 0.66 for T = 170 MeV and zero chemical potential, = 0. 

This is no longer the case when the baryon chemical potential is raised however. 
As an illustration, we trace the freeze-out line of Ref. ^Jj which reports freeze-out 
temperatures and chemical potentials for a wide range of collision energies y^. As we 
go to lower energies the temperature decreases slightly however the chemical potential 
rises sharply; as a result the population of baryons rises in comparison to the mesons. 
This leads to an increase in the coefficient Cbs as shown in the left panel of of Fig. ^ 
as a function of the chemical potential at the corresponding ^/s. 

The value of unity for Cbs in a quark-gluon plasma was obtained within the 
picture of a weakly coupled system. A more realistic value may be obtained from 
lattice simulations of QCD on the basis of off-diagonal susceptibilities ^21 ■ In these 
simulations, the partition function Z for QCD with three quark flavors is estimated. 
The input parameters are the temperature T, and the three chemical potentials for the 
the up, down and strange flavors (/i„, Ms)- Calculations are performed at vanishing 
chemical potentials. 

In this system the density of a flavor / is obtained as the derivatives of = 
\ogZ/V with respect to the chemical potential /j,/. The quark number susceptibilities 
are the second derivatives 
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It may be easily argued that the mean flavour densities are zero at vanishing 
chemical potential for all flavours: (u) = (d) = (s) =0. As a result, the coeSicient 
Cbs may be expressed in terms of the susceptibilities as 



, Xds Xus /-\ 

= 1 + -— . (5) 

Ass 



From Ref. [21 we obtain x^s/T^ = 0.53(1) and Xus + Xds = 0.00(3) at T = l.ST^. As 
a result we obtain C ~ 1 from lattice QCD, consistent with the estimate from a naive 
picture of a weakly interacting plasma. These results were obtained in a quenched 
approximation, but the effect of sea quarks is expected to be marginal above Tc |12| . 
We may thus surmise that the lattice system has Cbs ~ 1, which is consistent with a 
plasma of quasi-particle quarks. 

Recently there has appeared a model that purports to explain both the equation 
of state as obtained on the lattice as well as the large flow observed in heavy-ion 
collisions 0. The model describes the chromodynamic system as a gas of massive 
quarks, antiquarks, and gluons together with a myriad of their bound states generated 
by a screened Coulomb potential. In order to assess the consistency of this model with 
lattice calculations, we estimate the coefficient Cbs in such a scenario. 

Our estimates are based on Ref. |Sj and the assumption that the fluctuation 
coefficient Cbs is set at a temperature T = l.STc. In this model, the plasma 
at a temperature Tc < T < 3Tc possess massive quasi-particle excitations. The 
attraction between such colored states is modelled via a Coulomb potential with a 
Debye screening mass which is derived from lattice simulations. The temperature 
dependence of the quasiparticle and Debye masses is obtained by parameterizing data 
from the lattice. The attraction at such temperatures is still sufficient to produce a 
variety of bound states. For three ffavors of quarks there exist 749 such bound states. 
However, only the color-triplet sg and the color-singlet qs states (and their conjugates) 
are of relevance here. (The color-hexaplet sg states as well as the diquark states are 
very weakly bound and dissolve entirely at 1.5 Tc.) There are 2 7r-like (spin-singlet) 
and 6 p-like (spin-triplet) qs states, and 36 sg states. The abundances of these states 
are estimated in a grand canonical ensemble with vanishing chemical potentials. The 
qs multiplets carry no baryon number and thus contribute only to ag, hence drive Cbs 
towards zero, while more degenerate sg multiplet contributes also to ctbs and thus 
drives Cbs towards unity. The resulting value obtained at T = l.bTc, Cbs = 0.62, 
differs significantly from the value extracted directly from lattice QCD (see above), 
thus suggesting that the bound-state model, does not encode the same degrees of 
freedom that are pervasive in lattice simulations of QCD. 

Turning to more experimental considerations, we outline how such correlations 
may be actually measured in RHIC collisions. Such estimations are made with the aid 
of the Monte-Carlo event generator HIJING [Hj. The collisions occur in the center-of- 
mass frame at a -y/s = 2QQAGeV and we estimate the coefficient Cbs using the left 
hand side of Eq. . One measures all hadrons in a given event that lie within a range 
of rapidity from — lymaa:! < y accept < lUmaxl, and computes the total baryon number 
B, the total strangeness S and their products BS and S. One then calculates the 
event averages of all these quantities and resurrects the coefficient Cbs- These results 
are represented by the closed circles in the right panel of Fig. ^ 
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As the rapidity range of the acceptance is increased the coefficient remains more 
or less constant at a rather low value of 0.45. This is primarily due to the well 
know problem of baryon production in string fragmentation models. However as the 
rapidity range is increased we note a large rise in the coefficient. This is due to the 
fact that the net baryon number in the large rapidity regions is much larger than 
at mid-rapidity and as a result the number of baryons in the computation of the 
coefficient increases. This rise should be compared to the rise of the coefficient with 
baryon chemical potential in the left hand side of Fig. ^ Eventually, at very large 
acceptances, one begins to capture the entire system. Due to the strict conservation 
of baryon number and charge in such models, the fluctuations of BS drop rapidly as 
the acceptance encompasses the full system. There remains a residual fluctuation in 
the strangeness content due to weak decays. As a result Cbs ^ at large acceptance. 

In order to study the effect of the exact conservation of baryon number without 
the influence of the rather large net baryon number, we repeat the above study on 
JETSET (e"''e~annihilation to two jets) events ^3] where the net baryon number 
is zero. As would be expected, at the lower rapidity range of the acceptance, it is 
consistent with the results from HIJING, which also displays a very small net baryon 
number at mid-rapidity. The coefficient Cbs extracted from JETSET events shows a 
monotonous drop to zero as the acceptance is increased to encompass the full system. 
The value of Cbs for asymptotically small ranges of the acceptance may estimated 
as simply the ratio of the probability to observe a strange baryon to that of strange 
meson. 




Baryon chemical potential, Hg (MeV) Maximum rapidity accepted, y 



Figure 1. Calculations for the coefHcient Cbs = (SBSS) / {SS'^) in various 
scenarios. The dot-dashed line represents Cbs oil the lattice and in a quasiparticle 
plasma; the dashed line is the estimate for a hadron gas and for a Shuryak-Zahed 
(SZ) plasma. The filled circles represent the estimate from JETSET at center-of- 
mass energy = 20 GeV for various ranges of the maximum absolute rapidity of 
acceptance. The green triangles represent the results from HIJING simulations 
for Au-Au at 200AGeV for various ranges of acceptance. 

In these proceedings, we have proposed an experimentally observed ratio Cbs ~ 
—3{SBdS)/{6S'^), which may be used not merely to detect the prior existence of a 
chromodynamic phase of QCD, but indeed to decipher the strangeness and baryon 
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number carrying degrees of freedom prevalent in such an environment. Estimates 
based on a quasi-particle picture of such degrees of freedom tend to favor a value 
of unity for the coefficient. This is also consistent with the estimates from lattice 
simulations. In a more strongly interacting scenario, which possesses light bound 
states of quarks and antiquarks the coefficient is reduced to approximately 2/3 which is 
consistent with the coefficient from a hadronic resonance gas with no prior existence of 
a deconfined phase. Estimates from Monte Carlo models of string fragmentation place 
the coefficient at an even lower value of 0.45 nearly independent of the rapidity range 
of the acceptance up to ±2 units of rapidity. Thus measurements of the coefficient 
Cbs allow for a resolution of the baryonic and strange degrees of freedom in a quark 
gluon plasma. For further details see Ref. |14| . 

3. Strangeness at the SPS: /i > T 

In this section, the focus is turned to systems produced at the lower beam energies of 
the SPS. Deriving a freeze-out temperature T and a chemical potential fi from fits to 
particle ratios leads to a T ~ 165 — 160 MeV and hb — 240 — 300 MeV for a collision 
energy ^/s = 8 — 17 GeV. Most models @] indicate that at such large baryon chemical 
potentials the system may indeed undergo a first order phase transition from a partonic 
phase to a hadronic phase. As the collision energy is increased, the chemical potential 
tends to drop. As a result, the line of first-order phase transitions should eventually 
terminate in a second-order end point followed by the regime of continuous cross over 
as expected from lattice simulations at vanishing chemical potential. Current lattice 
estimates do indicate the presence of such a critical end point |3] . The location of this 
point is however not well established. 

The proximity between the expected first-order phase transition line and the 
chemical freeze-out line at these collision energies suggests that a quark-gluon plasma 
could possibly be formed at these energies. This system will then traverse the phase 
transition line during its subsequent evolution into a hadronic gas. In this section, 
we outline the observational consequences of a first order phase transition on the 
multiplicity and fluctuations of strange particles produced in such collisions. 

A universal feature of first-order phase transitions is the occurrence of spinodal 
decomposition, which occurs as a result of the presence of a a convex anomaly in the 
associated thermodynamic potential. This phenomenon occurs when bulk matter, by 
a sudden expansion or cooling, is brought into the region of phase coexistence. Since 
such a configuration is thcrmodynamically unfavorable the uniform system prefers to 
reorganize itself into spatially separate single-phase domains. These blobs of plasma 
immersed in an environment of hadrons tend to have a characteristic scale |15| . 
The essential feature of such a scenario is that if the breakup is sufficiently rapid, 
then whatever strangeness happens to reside within the region of the plasma that 
forms a given blob will effectively become trapped and, consequently, the resulting 
statistical hadronization of the blob will be subject to a corresponding canonical 
constraint on the strangeness. This will lead to an enhancement of the multiplicity 
of strangeness-carrying hadrons, as compared to the conventional scenario where 
global chemical equilibrium is assumed to be maintained through freeze-out. (This 
occurrence of an enhancement is qualitatively easy to understand, since the presence of 
a finite strangeness in the hadronizing blob enforces the production of a corresponding 
minimum number of strange hadrons.) The fluctuations in the multiplicity of strange 
hadrons, such as kaons, are enhanced even more, thus offering a possible means for 
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the experimental exploration of the phenomenon. Indications of an enhancement in 
the K-to-TT ratio has already been seen experimentally around 20 — 30 AGeV beam 
energy Furthermore, at the same energy, the fluctuation of this ratio are strongly 
enhanced if expressed relative to mixed events. 

Strangeness, as opposed to other conserved quantities such as charge and baryon 
number, is not brought in by the initial states in a heavy-ion collision. It is produced 
in the dense systems created in such collisions and thus the net strangeness in the 
full participant region is necessarily vanishing. There is however no constraint on 
the number of strangeness carriers e.g., the strange quarks and antiquarks or on 
the resulting strange mesons and baryons. For such an enhancement to occur the 
phase transition in question must necessarily be from a state with a greater number of 
strangeness carrying degrees of freedom (compared to the dominant carriers of entropy 
e.g., in a quark-gluon plasma) to a state with a fewer number of strangeness carriers 
as compared to the number of carriers of entropy {e.g., a hadron gas). The rapid 
expansion curtails the global chemical re-equilibration of the hadronic strangeness 
carriers emanating from different plasma blobs and results in a persistence of the 
overall enhancement of the mean number of strangeness carriers in the rarer hadronic 
phase. 

As in most statistical estimations of the multiplicities of produced particles, we 
consider a volume V in which the produced strange quarks and antiquarks are assumed 
to reach full global chemical equilibrium (at plasma temperature Tq, and baryon 
and charge chemical potentials ^b, I^q) under the canonical constraint of net total 
strangeness Stot = 0, i.e., we assume that the expansion is not rapid enough to offset 
the chemical or kinetic equilibrium of the produced quarks, antiquarks and gluons. 
The multiplicity distribution of the total number of strange quarks (identical to the 
distribution of anti-quarks) is given by the expression 



Where gg is the spin-color degeneracy of a strange quark and K^ix) = x'^K2{x) [K2{x) 
is the modified Bessel function of the second kind]. We denote the denominator of 
Eq. lO which represents the canonical partition function of a system of strange quarks 
and antiquarks with net strangeness zero as Zq''~ ^ = Zq. 

Following the dynamical scenario outlined above, we imagine that the system 
decomposes into p boxes of equal size Vq = V/p (here and in the rest of this section 
we will use blobs and boxes interchangeably). The number of ways of distributing N 
non- identical quarks into p boxes with ni quarks in box 1, n2 in box 2 etc. , is given 
by the combinatorial factor. 




The single quark (anti-quark) partition function ^ is given as 




(7) 




(8) 
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In the above equation, {n} refers to the set of occupation numbers of the various boxes 
or blobs i.e., rii, 77.2, ...Up. 

With the aid of the above expression and Eq. © we obtain the probability 
distribution of events with quarks and anti-quarks distributed in the p compartments 
according to the occupation vectors of quarks and antiquarks i.e., n'^ and n', as 



2!- 



(9) 



These blobs of plasma now undergo expansion and hadronization. The final 
hadronic populations achieve chemical equilibration within a volume Vh — x^q (^t a 
temperature Th) in the vicinity of the now well separated plasma blobs. The expansion 
factor is estimated to be x ~ 1-5 — 4. As a result, the net strangeness of a blob i 
{si = fii — Ui) is maintained post hadronization by the hadronic resonance gas. Thus 
a computation of the properties of the final hadronic gas requires the probability 
distribution of the net strangeness in the various blobs, 



' - E Pin'^,n'')l[6is, + n.-n.). (10) 

^ ni....np,ni,...np % 

In an experimental setup one never measures all the particles produced in the full 
phase space, but instead views a sub-sample of a given event. As a result, the focus 
will lie on the set of partially summed quantities. 



Ql{s,,...su)^\{zf^ ^ X{z%5{y.s,-s\ . (11) 

i=l Sfc + i,...Sp i \^~^ / 

Where Zf; represents the canonical partition function of only strange quarks and 
antiquarks in a single blob i with net strangeness Si. The above equation represents the 
probability distribution, that a system of strange quarks and antiquarks is distributed 
into p separate compartments and the strangeness in the first k compartments is 
specified. The strangeness contents of the remaining compartments is summed over 
all allowed values. 

Given the strangeness Sq in a compartment, the distribution of strange hadrons 
in that compartment may be computed using the canonical partition function for 
strange hadrons including the six classes of such hadrons with strangeness ±1, ±2, ±3 
and with total net strangeness Sq, 



^fr^'- n e&(e^-^-^o). (12) 

The effective single-particle partition function for the class of strange hadrons with 
strangeness S, C5, is a sum of the single-particle partition functions off all the hadrons 
belonging to this class. The general criteria for the inclusion a given strange hadron 
in the above equation will be that its mass be below 1680 MeV and its decay width 
be below 200 MeV. 
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The net strangeness in Eq. H12|l is given by the probability distribution 
Qp{So)/Zq computed above. The chemical potentials involved in the single-particle 
partition functions of the hadrons (Cfc) are in principle also influenced by the value of 
5*0 in the given blob. One may also note that the non-strange sector will be influenced 
by the traversal of the system through the first order phase transition line. In these 
proceedings we will ignore such complications and refer the reader to Ref. ^Jj for 
details. In what follows, we will assume the chemical potentials to be independent of 
the net strangeness in the blob. It will also be assumed that the non-strange hadron 
multiplicities and their fluctuations may be obtained by grand canonical estimates 
based on the measured freeze-out temperature and chemical potentials. 

As a result of the above simpliflcations, we may obtain the mean and variance of 
particles of strangeness k from a compartment of net strangeness Sq as JHl 

OQ OO 

The mean number of such particles from a single blob produced in a scenario as 
outlined above may be obtained by simply taking the convolution of the above equation 
with the probability distribution of strangeness, 

K) = ^E|ISQp(^o). (14) 

^" So '^So 

As all the different blobs or compartments are considered equivalent, we present 
results for the multiplicities and variances from a single such blob of volume 
Vq = 50 inr^ in the plasma phase in Fig. |21 For comparison, two other scenarios are 
presented. Besides the mean and variance that result from a spinodal decomposition 
during the phase transition, we also present the results from a simple grand canonical 
reference scenario. Using the baryon chemical potential as a control parameter, 
we obtain the freeze-out temperature Th from the fit to the data obtained in Ref. ^Hl ■ 
Subsequently, we perform a grand-canonical iteration to determine those values of /ig 
and /Lis that ensure Q = aB and S = 0, where a = 0.4 which is representative of Z/A 
for gold. The last is referred to as the restricted canonical scenario where the blob 
strangeness is always fixed to zero. 

The general trend in all the three scenarios is similar. As /.is is increased, the 

yield initially increases to balance the strangeness of the increasing number 
of strange baryons. As the freeze-out temperature begins to decrease noticeably, 
the hadron production generally decreases, leading to a decreasing behavior of the 

curve. The multiplicities are insensitive to the value adopted for the plasma 
temperature Tq, which governs the fluctuations in the blob strangeness Sq. For 
the variances, the overall behavior is qualitatively similar to the behavior of the 
averages, however, with some differences. First, in the restricted scenario (where 
only S'o = is included) the suppression of the variance is significantly larger. This is 
the well know effect of canonical suppression. Furthermore, at the larger values of ^b, 
where the net baryon density becomes significant, the grand-canonical variance suffers 
more from the decreasing temperature than the canonical variance. This important 
divergence is a result of the fact that the larger average baryon number implies a 
correspondingly larger variance in the baryon number and therefore also a larger 
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Figure 2. The average multiplicities and the A'+ variances as functions of 
the baryon chemical potential /ib in three three scenarios: 1) the standard grand- 
canonical treatment, 2) the canonical treatment in which the blob strangeness So 
is conserved through the hadronization (the mean multiplicities are not sensitive to 
the plasma temperature Tq), and 3) the restricted canonical treatment admitting 
only So = 0. 



variance in the strangeness. The enhancement is naturally reduced by reducing the 
plasma temperature as shown by the dashed line in Fig. |21 In both cases, the results 
from the scenario with a first-order phase transition with strangeness trapping exceeds 
those from the other scenarios. This is especially true for the variance. 
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